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Abst ract 


This  paper  presents  the  results  of  an  analysis  of  the  vibra- 
tional spectrum  of  the  solvated  alkali  metal  cations  in  dimethyl 
sulphoxide  (DMSO) . The  purpose  of  this  study  is  to  determine  more 
about  the  structure  of  the  primary  solvation  shell  and  to  formulate 
a clearer  picture  of  the  vibrational  states  of  these  species.  The 
ionic  systems  are  investigated  with  the  use  of  a potential  energy 
function  which  consists  basically  of  an  ionic  term  and  an  exponential 
repul s iui. . Dispersion  interactions,  where  they  apply,  are  accounted 
for  with  the  use  of  the  Morse  function.  We  show  that  for  systems  with 
cubic  symmetry  there  is  no  dependence  in  the  force  constant  upon 
electrostatic  terms.  Thus,  the  far  infrared  spectrum  provides  a 
probe  of  the  repulsion  interaction  which  operates  between  the  ion  and 
the  solvent.  On  the  other  hand,  the  breathing  modes  depend  upon 
all  contributions  in  the  potential  function.  Our  study  shows  that 
the  ions  of  lithium,  sodium,  and  potassium  can  be  regarded  as 
classical  point  charges.  Dispersion  effects  are  not  important  for 
these  ions.  The  far  infrared  spectra  of  the  ions  of  rubidium  and 
cesium,  on  the  other  hand,  depend  upon  dispersion  forces.  Finally, 
the  ions  of  lithium  and  sodium  mav  be  regarded  as  atomic  masses 
trapped  inside  an  infinitely  massive  solvent  cage.  For  the 
remaining  alkali  metal  ions,  reduced  mass  effects  are  important 
in  the  determination  of  the  vibrational  frequencies. 


1 nt  roduc  t ion 


A knowledge  of  the  mechanics  of  ionic  solvation  is  important 

for  the  understanding  of  a number  of  phenomena.  For  example,  the 

deformation  of  the  solvation  sheath  about  an  ion  contributes  to  the 

electrophoretic  effect . 1 An  adjustment  of  the  solvation  sheath 

about  an  ion  or  electroactive  species  often  contributes  to  the 

activation  of  an  electron  transfer  transition  in  an  oxidation- 

> 

reduction  reaction."  The  elastic  deformation  of  the  solvation 
structure  can  be  expressed  in  terms  of  t lie  vibrations  of  the  system. 
Thus,  a knowledge  of  the  vibrational  modes  of  solvated  ionic  and 
molecular  species  is  needed  in  order  to  understand  a number  of 
transport  and  kinetic  processes. 

Because  of  the  relative  simplicity  of  the  electronic  structure 
of  the  alkali  metals  and  their  cations,  they  are  potentially 
excellent  systems  to  investigate  in  order  better  to  understand 
many  of  the  features  of  the  electron  transfer  reaction. 

1 lect rochemi st rv  can  be  carried  out  on  the  alkali  metal  cations  in 
d i methyl  sulphoxide  (PMSO)  as  well  as  in  other  solvents.*  All  ol 
the  available  alkali  metal  cations  are  soluble  in  PMSO  to  some 
extent.  Thus,  for  these  systems  there  is  the  greatest  amount 
of  information  from  spectroscopic  experiments.  Thus, 
we  concentrate  our  attention  on  these  svstems  in  this  paper. 

In  particular,  we  are  interested  in  the  vibrational  spectra  of 
the  metal  cations  in  PMSO. 

Fdgell  et  al.**  and  Popov  and  his  colleagues^'’*'  have  carried 
out  extensive  investigations  of  the  far  infrared  vibrational 
spectra  of  the  alkali  metal  cations  in  PMSO.  Thus,  their  data  form 
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the  basis  ot'  our  investigation. 

4 So 

f.dgell  et  al.  and  Popov  ’ in  examining  their  data  have 
constructed  potential  energy  functions  in  order  to  try  to  determine 
the  force  constants  for  the  vibrations  in  these  systems.  Generally, 
their  analyses  have  been  based  upon  the  use  of  much  simplified 
models  of  the  primary  state  of  solvation.  In  particular,  both 
groups  have  examined  a linear  triatomic  model  in  which  the  ion 
is  sandwiched  between  two  solvent  molecules.  This  model  is  not 
adequate  in  that  it  does  not  usually  consider  enough  of  the  solvent 
which  surrounds  the  ion.  Specifically,  a linear  triatomic  model 
predicts  that  the  force  constant  will  depend  upon  an  electrostatic 
contribution  which  arises  from  the  ion-dipole  interactions.  In 
the  following  paragraphs  we  present  a refined  solvation  model  which 
can  account  for  more  complicated  states  of  solvation.  We  show  that 
whenever  the  geometry  of  solvation  satisties  simple  cubic  symmetry 
requirements,  there  will  be  no  electrostatic  contribution  to  the 
force  constant  for  the  far  infrared  vibrational  modes.  The 
conclusion  reached  by  Ihlgell  at  al.,  that  the  repulsion  interactions 
are  important  in  the  determination  of  the  magnitude  ot  the 
vibrational  force  constants  survives  in  our  treatment.  Indeed, 
it  is  strengthened.  We  assert,  in  fact,  that  the  far  infrared, 
dipole  allowed  transitions  of  the  ion  vibrating  within  the  solvent 
cage  provide  a direct  probe  of  these  repulsion  interactions. 

The  potential  energy  function  for  a system  of  an  ion  with  a 
shell  of  solvent  molecules  (as  dipoles')  can  be  constructed  from 
four  basic  parts.  The  first,  and  obvious,  part  consists  of  the 
electrostatic  interactions  which  operate  between  the  ion  and  the 
solvent  dipoles  which  surround  it.  The  second  part  consists  of 
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the  intcrmolecular  repulsion  interactions.  Third,  dispersion 
interactions  enter  to  a lesser  or  greater  degree  depending  upon 
the  degree  of  polarizability  of  the  ion  and  the  solvent.  Finally, 
if  the  ion  and  its  first  solvation  shell  are  immersed  in  a polar 
solvent  system,  it  is  necessary  to  consider  the  influence  ot  the 
remainder  of  the  solvent  system  on  the  central  ion. 

For  the  solvated  ionic  system  there  are  two  distinct  kinds 
of  force  constant  which  operate.  The  far  infrared,  vihrat iona 1 1 v 
active  transitions  depend  upon  the  oscillation  of  the  central  ion 
within  its  cage  of  solvent  molecules.  For  this  degree  of  freedom 
the  interaction  between  the  solvent  molecules  in  the  solvation  shell 
and  the  interaction  between  the  solvent  and  outer  solvent  molecules 
is  unimportant.  Indeed,  as  we  show  later,  electrostatic  interactions 
in  general  are  important  only  as  cubic  or  higher  order  corrections 
to  the  basic  spherical  harmonic  oscillator  potential  energy 
func t ion. 

On  the  other  hand,  the  symmetric  stretching,  Raman  active 
mode  depends  upon  a number  of  factors.  In  particular,  it  depends 
upon  the  elect rostat ic  terms  which  operate  within  the  primary 
solvation  system,  i.e.,  the  ion  and  its  first  solvation  shell, 
as  well  as  the  interaction  between  the  system  and  the  outer  solvent. 
The  secondary  interaction  between  the  primary  solvated  ion  and  the 
remaining  solvent  system  in  part  can  he  handled  in  terms  of 
electrical  image  interactions.  Specifically,  the  solvated  ion 
resides  in  a spherical  cavity.  The  charged  ion  and  the  solvent 
dipoles  all  have  images  in  the  surrounding  continuum  dielectric. 

In  addition,  there  is  a repulsion  force  which  acts  to  keep  the 
solvent  molecules  within  the  cavity.  As  a consequence,  although 
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the  breathing  inode  is  simple  with  respect  to  symmetry,  the  mechanical 
interactions  which  determine  the  strength  of  the  oscillation  are 
much  more  complicated  than  is  the  case  with  the  infrared  active 
vibration  of  the  central  ion. 

There  are  indeed  other  modes  of  vibration  which  are  possible. 

Many  of  these  modes,  however,  are  expected  to  be  of  very  weak 

intensity  and  low  frequency.  Nevertheless,  they  can  be  important  in 

the  activation  process  for  an  electron  transfer  reaction.  In 

particular,  the  simultaneous  stretching  along  one  axis  and 

the  contraction  along  another,  can  be  important.  This  type  of 

elliptic  deformation  has  been  discussed  before  in  connection  with  the 

7 

hydrodynamic,  collective  model  of  ionic  solvation.  It  represents 
effectively  a volume  conserving,  surface  oscillation.  In  this 
paper,  however,  we  do  not  examine  this  particular  vibrational  mode. 

On  the  one  hand,  there  is  no  direct  experimental  evidence  which  can 
be  attributed  to  such  vibrations.  On  the  other  hand,  it  may  well 
be  the  case  that  such  vibrational  modes  represent  collective 
oscillations  of  t he  secondary  solvation  shell.  In  spite  of  a 
lack  of  direct  experimental  evidence  for  the  existence  of  these 
vibrations,  they  are  nevertheless  important  to  the  electron  transfer 
process.  They  will  be  considered  more  in  detail  in  a separate 
paper . 

The  basic  form  of  the  interaction  between  the  ion  and  its 
collection  of  solvent  molecules  is  through  the  ion-dipole  inter- 
action. However,  in  order  to  be  able  to  describe  the  small 


vibrational  excursions  which  the  ion  can  make  away  from  its  equilibrium 
position  at  the  centre  of  the  solvation  cage,  it  is  necessary  to 
consider  the  expansion  of  the  interaction  in  appropriate  terms. 
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This  can  ho  done  in  terms  of  the  Car  1 son  - Rushhrooke  formula. 

The  reason  for  using  the  Car  1 son - Rushhrooke  expansion  is  the 
following.  In  a system  in  which  the  ion  moves  with  respect  to 
a fixed  set  of  solvent  dipoles,  it  is  necessary  to  he  able  to 
specify  the  interaction  between  the  ion  and  each  of  the  dipoles 
with  reference  to  one  point.  Clearly,  for  an  ion  with  some  torm 
of  symmetric  solvation,  the  reference  point  for  the  ion  is  at  the 
centre  of  symmetry.  The  use  of  the  Carlson  Rushhrooke  expansion 
formula  [see  appendix  l for  details]  enables  us  to  express  the 
ion-dipole  interaction  in  terms  of  coordinates  which  are  appropriate 
to  the  analysis  of  the  vibrations  of  the  ion  in  the  cage.  This 
has  not  been  done  in  other  treatments.  Our  use  of  this  type 

of  expansion  allows  us  to  consider  the  types  of  caged  ion  vibrations 
which  occur  in  condensed  systems,  but  not  necessarily  in  the 
gas  phase. 

In  the  following  section  we  present  the  various  forms  of  the 
potential  energy  function.  The  force  constants  which  are  derived 
from  these  functions  are  summarized.  The  details  of  many  of  the 
calculations  are  presented  later  in  this  paper  in  a series  of 
appendices.  There,  a number  of  technical  points  are  developed. 

These  matters,  pertaining  to  computation,  are  important  in  the 
analysis  of  the  activation  of  the  electron  transfer  reaction. 

The  application  of  the  force  constants  to  the  analysis  of  the  spectrum 
of  the  solvated  alkali  metal  cations  in  PMSO  is  considered  in  the 


t h i rd  sec  t i on . 


Potential  I ih'i  o i unc t ion 


and  \ ibrat  ion.il  force  I’onstants 


Iu  this  section  we  present  potential  energy  functions  for 
two  haste  types  of  solvated  ionic  system:  a four  - coord i nat ed  , 
tetrahedral  solvation,  and  a six  coordinated,  octahedral  solvation, 
'lent  ion  is  made  ot  t lie  e i ji  h t coo  rd  i na  t ed  , cubic  solvation  case. 

lor  an  arrangement  of  solvent  dipoles  tin  a sense  to  ensure 
an  attractive  interaction)  at  the  vertices  ot  a regular  tetra 
hedron,  t he  larlson  Kushbrooko'  expansion  yields 

\ t I 1 M j.  liW"  l V , : ( r ) V , , 1 r ) U r/Rl  ' ] ] ( 1 ) 

l\  * ^ 

in  which  V,  ir)  is  the  spherical  harmonic  function.  five  vector 
r = t.  x , v , : ) det’mes  the  position  of  the  ion  with  re  to  fence  to  the 
centre  of  symmetry.  I he  distance  R is  the  distance  between  the 
centre  of  symmetry  and  the  solvent  molecules.  In  the  regular 
array,  each  molecule  has  the  same  value  ot  K.  The  molecular  dipole 
moment  is  jjiven  bv  n. 

\n  octahedron  with  dipoles  at  the  vertices  yields  the  following 
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Another  form  which  is  possible,  but  probable  not  common , 
places  the  dipoles  at  the  vertices  of  a regular  cube.  In  this  case 
the  first  term  past  the  .•eroth  order  contribution  depends 
Olr'1 ) . 


upon 


- s - 

It  is  obvious  that  the  equilibrium  position  is  located 
at  r = 0,  i.e.,  at  (0,0,0).  It  is  equally  clear  that  in  both  eqn  (1) 
and  (2)  there  no  contribution  to  the  potential  of  order  0(r“). 
Hence,  when  the  harmonic  oscillator  limit  is  sought  (through  the 
use  of  a Taylor  series  expansion)  in  the  complete  potential  energy 
function,  there  will  he  no  electrostatic  contribution.  This  is 
not  true  if  the  solvation  is  linear.  Tor  linear  solvation,  however, 
two  distinct  caged  ion  vibrational  modes  are  predicted;  this  is 
considered  further  in  an  appendix.  The  non - i nvol vemcnt  of  electro- 
static terms  in  the  force  constants  anpl  ies  only  to  the  infrared, 
dipole  allowed  transitions  in  a point  charge  sistom. 

There  are  several  acceptable  expressions  for  the  usual lv  short 
ranged,  repulsion  interaction  which  operates  between  atomic  and 

- 1 2 

molecular  species.  Although  the  repulsion  term  (R  -)  in  the 

l.ennard -Tones  potential  is  perhaps  best  known  and  most  widely  used, 

n 

an  exponential  form,  used  extensively  by  Horn,  can  also  be 
employed.  There  are  obvious  computational  advantages  to  the 
use  of  the  l.ennard -.Jones  potential  for  many  applications.  However, 
for  the  cent  ro  - symmet  r i c systems  we  investigate  here,  the  l.ennard  - 
• Jones  function  lias  the  disadvantage  of  not  allowing  an  appropriate 
I.ap  1 ace  * type  expansion.  Thus,  for  the  repulsion  which  operates 
between  species,  we  make  use  of  the  exponential  representation. 

The  dispersion  interaction,  as  it  is  contained  in  the  Morse 
potential  function,  is  also  given  in  an  exponential  form. 

The  Born  repulsion  has  the  simple  exponential  form 
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wlu'iv  H :i  iu!  P .ill'  constants.  When  dispersion  is  a I so  considered 
the  Ho  rn  t o rm  is  expanded  into  the*  Morse  funct  ion;1*1 


I exp  [ - C r - r ) ] | ? 


wliori',  a s usual,  n is  the  il  i s soe  i a t i on  energy  ami  k is  a const  a. 
with  the  units  of  inverse  ili  stance. 

I he  lap  1 ace  expansion  of  the  exponential  function  (o')  is 
j;  i ven  hv  (see  appemlix  3) 


Vr  1,,H  ‘’*Pt»-e/t’)(i>;jfmCr)vCin(,n{(2l+niflr/p)kt(R/p) 

♦ (r/(Oie  + 1 (r/oHe(R/pl  (R/p)i(lr/(')k{  + 1(R/p)}.  (5) 

I lie  I unctions  i ( \ ) ami  k^(x)  are  respectively  the  modified 
spherical  Hessel  functions  of  the  first  and  second  kinds.11  This 
i xp. ins  ion  is  val  id  on  I \'  as  long  as  r < It.  A more  general  expression 
1 an  easily  and  obviously  he  written  (see  appemlix  3). 

for  a lour  coordinated,  tetrahedral  system  the  contribution 
I rom  ei|ii  (31  through  (3)  to  the  system  potential  energy  function 
is  given  hv 
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Hi' t «*il  ahovo,  a I f lioutth  tho  lort’o  Oi’itstant  ili'os  ni't  tlopoiitl  upon 

anv  o I oo  t i os  t a t 1 1 oon t i i hut  i oils  , thoro  avo  contributions  from 


the  repulsion  terms  ami  from  the  dispersion  interaction  (as  it  is 
expressed  in  the  Morse  function). 

In  the  limit  in  which  dispersion  interactions  play  little  or 
no  role,  and  where  there  are  only  repulsion  contributions,  the 
expressions  for  the  force  constants  in  the  harmonic  oscillator 
limit  are  the  following:  the  four  coordinated,  tetrahedral  system, 

k,  = 14H/.V*  HI  2e/Rtj;  (9a 

the  s i x coordinated,  octahedral  system, 

k()  = CU/p-’lU  :p/Re);  (9b 

and  the  eight  coordinated,  cubic  system, 

ks  = (8H/.VZ)U  - 2p/R  ) . (9c 

In  contrast,  a system  potential  which  is  a simple  combination 
of  an  electrostatic  term  and  a Morse  function  ( i . e . , dispersion 
included)  yields  the  following  expressions  for  the  force  constants 

V,  = 


8 iK  ■'/:■> 


(10 
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] 

1 1) 

Moelwyn-liughes  points  out  that  for  small  deviations  about 
an  equilibrium  position  the  henna rd -Jones  potential  and  the  Morse 
potential  are  equivalent.  Thus,  k can  be  replaced  bv  t>/Re,  and 
eqn  11)  can  be  expressed  as 

h,  = 9oI)/ir  Cl  la) 

4 e 

k = 144P/R2  (lib) 

b e 

The  force  constants  are  seen  to  vary  as  !!/  as  the  equilibrium 
solvation  radius  changes.  There  is  also  a radial  dependence  in 
the  dissociation  constant  P.  Altogether,  as  we  shall  show  in 

. g f ] 

the  next  section,  the  force  constants  depend  upon  as  Rp  . 

Finally,  it  is  necessary  to  consider  the  fact  that  the  solvated 
ion  occupies  space,  effectively  in  a cavity,  inside  the  bulk  oi 
a continuum  dielectric.  However,  because  we  have  found  that  the 
image  contributions  to  the  breathing  mode  force  constants  amount  only 
to  a small  change  in  the  frequencies,  we  defer  our  discussion  of 

f 1 

these  contributions  to  an  appropriate  appendix.  It  is  clear  that 
as  far  as  we  are  concerned  with  the  harmonic  oscillator  limit 
for  the  infrared  active,  ionic  vibrations  within  the  primary 
solvation  cage,  the  image  effects  (also  electrostatic)  will  be 
non - cont r i but i ng . The  image  contribution,  as  we  will  show,  depends 
importantly  upon  the  size  of  the  cavity  within  which  the  ion  resides. 

If  the  cavity  dimensions  change,  then  there  can  be  an  important 
image  cont r i but i on . [We  do  not  consider  this  effect  in  this  paper.) 

The  system  potential  energy  function  which  is  a combination 


of  the  electrostatic,  repulsive,  and  dispersion  terms  can  be 


r 
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differentiated  with  respect  to  R,  the  solvation  radius,  to  obtain 
the  force  constants  for  the  breathing  modes.  In  addition  to  the 
terms  which  have  been  mentioned,  the  breathing  modes  depend  upon 
the  solvent  d i po le-d ipole  interactions  as  well  as  upon  an  inter- 
molecular  repulsion  of  the  Horn  type.  For  a tetrahedral  system, 
we  can  write 

V.  = 4—  ♦ 4 b exp  [ - (R  R )/p]  + Ob’  exp  f - (2/77?R-R  • )/p  • ] 

’ 1 R • 1 

♦ (lS/4)/VTlvi?/R3)  (12) 

in  which  the  constants  b’  and  p'  account  for  the  possibility  that 
the  i nt ermo 1 ecu  1 a r interaction  may  differ  from  the  ion-solvent 
interaction.  We  assume  that  the  two  types  of  interaction  are 
roughly  similar.  Therefore,  the  force  constant  for  the  breathing 
mode  is  obtained  directly  (note,  this  term  does  iuvt_  include  any 
account  of  the  electrical  image  interactions) 

k,  = - 24  — + 45/5/7  111  + 20b/ p 2 . (13) 

bl  R1*  R s 

e e 

In  the  next  section  we  consider  the  application  of  these 
results  to  the  examination  of  t he  vibrational  spectra  of  the 


alkali  metal  cations  in  I'MSO. 
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Piscussion  of  tho  Solvation  Spectra 

The  data  of  Maxey  and  Popov''’  on  the  far  infrared  spectra  of 
the  alkali  metal  cations  in  PMSO  are  reproduced  in  Table  1.  These 
experimental  data  form  the  basis  of  our  discussion  . 

A preliminary  analysis  of  these  data  is  needed  in  order  to 
extract  the  force  constants.  We  assume  that  the  harmonic  oscillator 
limit  appl ies : v i z . , 

!/A  = ^^717  (14) 

is  the  relationship  between  the  frequency  (in  wave  numbers),  the 
mass  n and  the  force  constant  k;  ('  is  the  velocity  of  light. 

At  this  point  it  is  necessary  to  specify  the  form  of  the  mass 
quantity  which  is  to  be  used  in  eqn  (14).  In  general,  the 
expression  for  the  reduced  mass  of  the  ion  (mass  m.)  which 
oscillates  with  respect  to  e neighbouring  solvent  molecules 
(each  m ),  is 

u = cm  m . / ( c m + m . ) (15) 

sis  i 

In  the  limit  of  large  c,  the  reduced  mass  is  the  same  as  the  ionic 
mass.  This  limit  also  applies  to  the  case  of  an  infinitely  massive 
solvent  system.  In  other  words,  the  ion  vibrates  within  a cage  of 
rigidly  fixed  solvent  molecules. 

With  the  reduced  mass  as  specified  by  eqn  (15),  it  is  possible 
to  obtain  values  for  the  force  constants  for  the  ions  under  the 
assumption  of  various  degrees  of  coordination.  The  force  constants 
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aro  l i > t o*.l  in  Tahlo  -.  From  tin*  tablt*  it  appears  that  the  lithium 

cation  has  a frequency  which  is  at  odds  with  the  rest  for  any 

type  of  coordination.  Hie  cations  of  sodium  and  potassium  seem 

to  pair  and  the  cations  of  rubidium  and  cesium  likewise  seem  to 

form  a pair.  We  show  shortly,  however,  that  the  frequency  of  the 

lithium  system  is  wholly  consistent  with  the  frequencies  of 

sodium  and  potassium.  The  frequencies  of  the  rubidium  and  cesium 

systems  are  distinct  from  t he  rest.  In  any  case,  it  is  clear 

from  the  table  that  it  is  not  reasonable  to  assume  that  the  force 

constants  are  uniform  for  the  ionic  systems  of  the  group.  Such 

4 

an  assumption  was  made  by  Maxey  and  Popov. 

In  order  to  make  the  comparison  between  t he  data  and  the 

formulae  for  the  force  constants  presented  in  the  last  section,  it 

is  necessary  to  specify  the  values  ot  the  solvation  radii.  We 

assume,  as  have  Abraham  and  1 iszi,  “ that  the  solvation  radius 

can  be  expressed  as  the  sum  of  t lie  ionic,  Goldschmidt  radius 

l 2 

and  the  radius  of  the  solvent.  Abraham  and  I. iszi  “ calculated  the 
solvent  radius  from  the  liquid  partial  molar  volumes  with  the  use 
of  the  Steam  Tyring  formula.*'  The  values  of  the  radii  are 
listed  in  Table  A. 

first,  we  compare  the  experimental  results  with  t lie  theoretical 
expressions  for  the  force  constants  obtained  under  the  assumption 
that  dispersion  plays  a minor  role.  Thus,  the  force  constants 
of  cqn  (Pi  apply. 

When  the  simple  exponential  repulsion  is  used,  the  two  para 

1 4 

meters  R and  p must  be  specified.  The  work  of  Huggins  anil  Mayor 
on  the  alkali  halide  crystals  indicates  that  o for  those  systems 
is  essentially  a constant.  It  seems  reasonable  to  make  a similar 


w 
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assumption  hero  for  the  solvated  ionic  systems.  There  is  no  reason, 

however,  to  expect  that  it  should  he  the  same  constant. 

The  parameter  B can  he  eliminated  in  favour  of  other  quantities 

which  are  more  easily  determined,  vi:.,  the  ionic  charge,  the 

solvent  dipole  moment,  and  the  equilibrium  separation  between  the 

o 

ion  and  the  solvent  molecules,  following  Born's  treatment  ol  the 
crystalline  alkali  halides,  we  write 

\ » -cU'u/R*  - B exp(  - (R-R  )/pl  1 (lb) 

where  again  c is  the  coordination  number.  The  simple  equilibrium 
condition  dV/dR  - 0 yields 

B * Jeup/R^ . 117) 

- o 

We  find  that  an  optimal  value  of  p is  1.9  \ 10  ' cm.  With 
the  use  of  this  value  of  p together  with  the  values  of  the 
solvation  radii,  we  determined  the  values  of  the  force  constants 
which  are  listed  in  Table  4. 

The  results  shown  in  Table  4 indicate  a clear  grouping  of  effect 
Indeed,  the  theoretical  results  indicate  that  lithium  belongs 
with  sodium  and  potassium  in  its  behaviour.  The  force  constants 
and  the  associated  frequencies  are  reasonable  accurately  predicted 
for  the  series  of  ions  lithium,  sodium,  and  potassium  under  the 
assumption  of  four -coord inat i on . It  is  clear  from  the  results 
that  s i x -coord inat ion  for  these  three  ions  is  inconsistent  with  the 
experimental  data.  It  is  worth  noting  from  the  table  the  fact 


1 


t h .•  t w i t h t he  use  o i t hi*  infinite  m.i  s s assumption  tor  the  solvent, 
tlu*  best  agreement  w 1 1 h experiment  is  found  tor  lithium.  That  is, 
lithium,  which  has  .1  small  mass  aiul  verv  small  polarizability, 
shows  nee.  1 i g i b 1 e dispersion  am!  reduced  mass  effects.  At  the  other 
extreme,  potassium,  winch  lias  a mass  about  half  the  mass  of  the 
solvent,  begins  t o show  a definite  dependence  upon  the  choice  of 
the  appropriate  reduced  mass.  Ihe  results  for  sodium  indicate 
that  this  ion  is  better  handled  with  the  assumption  that  the  cage 
is  rigid.  Ihe  value  o f .'i'(>  cm  * for  the  vibrational  frequency 
compares  favourably  with  the  experimental  value. 

Turning  to  the  ions  of  rubidium  and  cesium,  it  seems  that 

tour  coordination  is  inconsistent  with  the  experimental  data. 

The  assumption  of  six  coordination  brings  the  theoretical  results 
more  into  an  agreement  with  the  experimental  results.  However, 
the  unmistakable  fact  that  the  t heoret i ca 1 1 v predicted  force 
constants  decrease  from  rubidium  to  cesium  is  at  odds  with  the 
experimental  fact  that  the  force  constants  increase  in  t hi  s 
direction.  The  fault  lies  with  the  neglect  of  the  polarizability 
of  these  larger  lonw  systems. 

\ s \,o  have  indicated  in  the  last  section,  the  Morse  potential 
can  be  used  to  bring  an  account  of  the  dispersion  interactions  into 
the  analysis.  Tor  small  amplitude  oscillations  it  is  possible 
t o argue  the  t runs  forma  t ion  from  the  l.cnnard  .tones  to  the  Morse 
potential.  The  quantity  P in  the  Morse  potential  is  the  dissocia- 
tion energy  for  a single  solvent  ton  couple.  hv  means  of  an  argument 
outlined  bv  Moelwvn  Hughes,'1  it  is  possible  t o identify  the 


dissociation  energx  further  in  terms  of  the  1 oinlon  force  which  operates 
between  the  ion  and  the  solvent  molecule.  Thus,  the  dissociation 


- I s - 


energy  can  be  expressed  as 


V • -i—  118) 

:r* 

e 

where  C is  a constant  which  can  be  identified  by  comparing  eqn  (18) 
with  the  London  interaction:^*'* 


31  1 it  a . 

v.  , = ■*-  ^ -1 . UiO 

London  ,R* . , 

si 

In  t li  i s expression  1 and  1.  are  respectively  the  ionization 
potentials  of  the  outermost  electron  in  the  solvent  and  the  ion. 

The  polarizabilities  of  the  solvent  and  the  ion  are  expressed  by 
the  terms  as  and  a..  In  eqn  (19)  it  is  apparent  that  the  ioniza- 
tion potential  for  the  outermost  electron  in  the  ion  is  a quantity 
which  general lv  will  be  much  larger  than  the  ionization  potential 
for  the  outermost  electron  in  the  solvent.  Therefore,  we  can 
simplify  eqn  (19)  by  writing 


V 


l.omlo  n 


, l a a - 
O SSI 


CO) 


The  dissociation  energy  now  is  expressed  as 


Ii 


T 


l a a 
s s 


i 


R* 


e 


Cl) 


When  this  expression  for  0 is  substituted  into  the  expressions  for 
the  force  constants,  we  find 


k 


4 


I a a 
s s 


i 


R" 
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and 


(23) 

is  difficult  to  determine  the  ionization 
solvent  molecule.  Information  about 
the  ionization  potential  of  DMSO  is  scarce.  Nevertheless, 
these  expressions  for  the  force  constant  still  can  be  used  in  a 
slightly  modified  sense.  In  particular,  we  are  interested  to  see 
if  it  is  possible  to  predict  the  value  of  the  force  constant  for 
another  species  (in  a homologous  series)  given  the  experimental 
value  of  an  ion  in  the  series.  To  do  this,  we  simply  form  the 
ratio  of  the  force  constants: 


1 art- 

k = 108  ■ — 1 

b P « 


In  some  instances  it 
potential  of  a particular 


k 2 / k i - (Ri/Ril^fai/ui) 


(24) 


This  ratio  holds  as  long  as  both  species  have  the  same  coordination. 
The  comparison  between  force  constants  for  four-  and  s ix -coord i nat ed 
spec i e s is  the  foil ow i ng : 


k 


(4) 


l‘\tOA\4)K 


(25) 


In  order  to  test  eqn  (24)  we  make  use  of  the  ionic  refractions 
(related  to  the  polarizabilities  by  a common  constant  factor) 
given  by  Smyth;*1'  they  are  reproduced  in  Table  5.  The  results  are 
shown  in  Table  b.  It  can  be  seen  from  the  Table  b that  for  sodium 
and  potassium,  as  determined  from  lithium,  there  is  no  agreement 
with  experiment.  The  inclusion  of  any  account  of  dispersion 


On  t he  other 


JO 

interaction  is  not  appropriate  tor  these  systems, 
hand , the  agreement  between  the  calculated  and  experimental  values 
tor  t tie  force  constants  tor  cesium  as  calculated  from  the  rubidium 
force  constants  in  vorv  good.  The  increase  in  the  value  of  the 
mice  constant  for  all  values  of  the  coordination  is  correctly 
p red  1 c t ed  . 

! qn  (J4)  reasonably  accurately  predicts  the  force  constant 
for  rubidium  from  potassium;  this  implies  the  coordination  remains 
the  same.  1 he  ratio  ot  the  force  constants  (9b)  to  (9a)  is  1.5. 

\s  we  discussed  earlier,  this  ratio  also  gives  a reasonably  accurate 
prediction  of  the  rubidium  force  constant  from  that  of  potassium. 

In  this  case,  however,  the  agreement  is  not  quite  as  good  as  that 
found  with  eqn  (Jl).  On  the  other  hand,  if  we  use  eqn  (J5)  to 
attempt  to  predict  the  change  in  force  constant  in  going  from 
potassium  to  rubidium,  the  predicted  force  constant  is  30,000  dvnes/cm 
too  large.  As  Svmth  notes/’  the  ionic  refractions  are  not 
very  accurate  quantities.  Thus,  there  may  be  some  uncertainty 
associated  with  their  use  in  tins  case.  In  any  event,  the 
coordination  of  the  ions  of  rubidium  and  cesium  remains  a cloudy 
issue. 

We  summarize  the  far  infrared  vibrational  analysis  as  follows. 

1 ithium  and  sodium  behave  as  classical  point  charges  which  vibrate 
within  rigid  solvent  cages.  The  magnitudes  of  the  force  constants 
for  these  two  systems  are  determined  entirely  In-  the  form  of  the 
short -ranged  repulsion  interactions  which  operate  between  the  ion 
and  the  solvent  molecules.  Potassium  also  behaves  as  a classical 
point  charge.  Its  force  constant  depends  only  upon  the  short 
ranged  repulsion  interactions.  However,  potassium  has  a mass 


*i«IW 


:i 

which  is  approximately  halt  tho  mass  ot  the  solvent.  lienee,  it 
is  appropriate  to  consider  the  reduced  mass  in  determining  the 
value  of  the  ionic  vibrational  frequency.  Rubidium  and  cesium 
appear  to  he  s i x -coord  i nat  oil  in  PM  S 0 , but,  as  noted  immediately 
above,  this  cannot  be  a definitive  coordination  on  the  basis  ot 
the  spec t roscop i c data  alone.  These  ions  of  rubidium  and  cesium 
depend  important ! v upon  the  oporat ion  ot  the  dispersion  forces  ot 
the  l.ondon  tvpe.  I'he  fact  that  the  force  constant  increases  in 
going  from  rubidium  to  cesium  is  accounted  tor  accurately  when  the 
polarization  of  the  ions  is  taken  into  account.  In  effect,  the 
importance  ot  these  dispersion  interactions  simple  points  to  the 
fact  that  these  ions  cannot  be  reyarded  as  simple  classical  point 
charges.  Quantum  mechanical  effects  enter  and  must  be  considered. 

I mat lv,  we  turn  to  a cons i derat i on  ot  the  breathing  modes 
for  these  svstems.  Mnfort unate  1 v , the  experimental  picture  is  not 
clear.  1 il'lOj  lias  been  examined  in  acetone.  In  the  Raman 

spectrum  there  is  a shoulder  on  the  act ive  1 ine  ot  acetone  which 
is  about  14  cm  ' greater  than  the  acetone  line.  It  we  assume 

that  this  is  due  to  a symmetrical  stretching  breathing  mode,  then 
for  acetone  we  calculate  the  force  constant  to  be  . 2 x 10^  dynes/cm 
Using  the  general  formula  for  calculating  the  t requency  tor  a 
symmetric  stretching  mode , ^ 


v7  = k/m  (,2(il 

s 

where  m.  is  the  mass  of  the  solvent  molecule  in  a tetrahedral  XY^ 
svstem , we  determine  that  the  frequency  should  be  134  cm 
Agreement  with  the  assumption  that  the  14  cm 


shoulder  is  a 


symmetric  stretch  is  apparently  not  pood . However,  it  is  not 
necessarily  reasonable  to  expect  to  be  able  to  use  the  simple 
solvent  mass  in  this  case.  The  symmetric  stretch  is  a volume 
non  conserving  mode.  Thus,  the  expansion  against  the  remainder  of 
the  system  may,  in  effect,  amount  to  an  expansion  of  a much  larger 
e f t e c t i v e mass.  \n  effective  mass  of  about  5384  gm  gives  agreement 
assuming  that  the  value  of  the  force  constant  is  accurate.  In  any 
case,  the  experimental  situation  is  not  definitive  enough  for  us 
to  make  any  clear  statement  about  these  breathing  modes. 

One  conclusion  we  can  make,  however,  is  that  in  the  image 
approximat ion , the  effect  of  the  image  terms  on  the  breathing  modes 
is  insigni f i cant . Hie  variation  of  the  force  constant  versus  the 
ionic  solvation  radius  is  shown  for  the  cavitv  and  the  cavity 
plus  its  image  cent r i but  ions  in  figure  1.  \s  can  be  seen,  there 
is  very  little  difference.  Consequent  1 y , we  believe  that  these 
terms,  which  are  an  additional  burden  to  calculate,  can  be  ignored. 
It  should  be  noted , nevertheless,  that  in  certain  applications, 
in  particular  in  the  calculation  of  electron  transfer  activation 
energies,  these  image  contributions  need  to  be  considered.  For 
the  vibrational  problem  alone,  they  are  small. 


Ibis  work  was  supported  bv  the  H.  S.  Office  of  Naval  Research, 


\rlington,  Virginia,  USA. 


-23- 

Ap pend  Lx  1 : Electrostatics  of  ionic  solvation 


The  derivation  of  the  Carlson-Rushbrooke  expansion  has  been 

8 19 

discussed  in  detail  by  these  authors.  Ruehler  and  Hirschfelder 

were  the  first  to  point  out,  however,  the  fact  that  the  original 

expansion  is  not  a general  formula.  Further  discussion  and 

generalization  of  the  original  expansion,  therefore,  has  been 

19  20  21 

given  by  Ruehler  and  Hirschfelder,  * by  Ruedenberg,  , by  Sack, 

2 2 

and  by  Kay,  Todd,  and  Silverstone.  For  our  purposes,  the  original 
expansion  suffices.  We  are  therefore  limited  to  a region  of 
applicability  such  that  the  sum  of  the  distance  magnitude  of  the 
dipole  (viz.,  u/e)  and  the  distance  of  the  ion  from  its  equilibrium 
position  is  less  than  the  distance  between  any  solvent  molecule  and 
the  centre  of  symmetry.  This  limit  is  made  mathematically  clear 
below . 

The  general  formula  is  given  by 


1/r 1 2 = l ( - 1 ) 
l ,m 

l ' ,m ' 


a a ' 

r i r 2 
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d+n ' +l 
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and  R >ri  + r2.  The  R-coefficient  is  given  by 
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The  coo  ill  i na  t e svstom  applicable  to  eqn  (1.1)  is  shown  in  liy.  J. 

W'hon  one  of  the  interacting  particles  is  a point  dipole  of 
magnitude  vi  = eis  , for  example,  the  expression  can  he  applied  t o 
the  types  of  systems  we  have  discussed  in  the  main  hodv  of  this  paper. 
I'he  remaining  coordinate  iq  now  specifies  the  location  of  tlu* 
ion  with  respect  to  the  or i yin  as  shown  in  1 iy.  J.  This  is  precisely 
the  type  of  coordinate  which  is  needed  in  order  to  specify  the 
vibrational  excursions  of  an  ion  awav  from  its  centre  of  symmetry 
in  the  solvated  system. 


for  a system  \>  i t h tetrahedral  solvation,  we  outline  the  steps 
necessary  to  determine  the  form  of  the  electrostatic  contribution 
to  the  system  potential  eneryv  function. 

The  It  coefficient  in  this  case  is  specifically 
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(1.4) 


The  b coefficients  are  listed  in  Table 

To r the  tetrahedral  system  it  is  convenient  to  place  the  solvent 
dipoles  at  the  following  locations: 

(1,1,1)  (1,1,1)  (1,1,1)  (1,1,1) 

Ity  means  of  scaliny,  t lie  solvation  radius  can  be  inserted  into 
the  potential  eneryv  function.  In  the  same  order  as  locations  are 
listed,  the  direction  cosines  for  the  dipoles  are  chosen  to  be 
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tlu'  I'omauuiu',  toruis  tollou  oasilv  hv  svwrnot  rv  : 


Tlu'  not  soi  oml  oiiloi  i out  r i blit  i on  is  t ho  sum  ol  (l.S)  to  (I, SI. 

Iliiis  , thoro  is  no  soooml  onlor  oon t r i but  i on . 

I'ho  no\t  oiiloi  i out  i i but  i on  is  t hi  r il  onlor  in  r . I'ho  rosult 
is  i;  i roil  in  tlu*  to.vt  ol  t ho  papor,  oijn  (I).  I'ho  valuos  ol  him'.iii) 
nooiloil  in  that  oaloulntion  aro  list  oil  m I’ablo  S. 

i'ub  to  solvation  (oio.lit  ooonl  i na  t ion)  oan  ho  shown  oasilv  to 
violil  a /on'  tli  i nl  onlor  to  nil.  I'ho  vortioos  ol  a rot’ll  I nr  oubo 


can  In*  rega  rded  in  terms  of  the  vertices  of  two  i liter -penet  rat  i ng , 
identical  tetrahedra.  The  rotation  about  any  t'^  axis  takes  one 
tetrahedron  into  the  other.  The  sum  of  the  contributions  from 
the  two  tetrahedra  cancel.  Thus,  the  first  order  contribution 
in  the  cubic  system  is  of  order  0(H). 

The  remaining  general  type  of  solvation  to  consider  is  that 
of  the  linear  "triatomic"  type  considered  before  by  bdgell,  et  al  . , 
and  bv  Maxev  and  Popov."’  Although  the  authors  indicate  that  the 
use  of  tins  model  is  to  an  extent  suggestive  only  of  the  types 
of  forces  which  need  to  be  considered,  Maxev  and  Popov'"’  give  linear 
solvation  considerable  credence  on  the  basis  ot  their  NMR  solvation 
s t ltd  v . Their  analysis  suggests  that  in  PMSO  l.i+is  stronglv 

associated  with  two  molecules  of  the  solvent.  Their  analysis  also 
suggests  that  for  the  ions  of  sodium  through  cesium,  the  coordination 
is  less  we  11  deli ned . 

There  are  difficulties  associated  with  the  use  of  a linear 
solvation  model  which  have  to  do  with  the  problem  of  the  agreement 
between  theory  and  experiment.  In  particular,  the  linear  solvation 
model  yields  two  non  vanishing  cont  r i but  ions  of  order  Oil**")  to 
the  system  potential  energy  function.  One  ot  these  contributions 
is  an  exponential  repulsion  term,  and  the  other  is  electrostatic. 
Specifically,  the  electrostatic  contribution  can  be  written  as 
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omitting  the  zeroth  order  contribution.  The  repulsion  contribution 
has  a similar  functional  dependence  on  the  coordinate  variables. 


It  is  clear  from  t he  Conn  of  this  potential  function  that  it  already 
contains  t lie  appropriate  harmonic  oscillator  potential  terms. 

Moreover,  it  is  clear  that  there  are  two  infrared  active  vibrational 
modes  for  the  solvation  system.  The  vibration  atom;  the  axis 
of  symmetry  o one  mode,  and  the  motion  of  the  ion  perpendicular 
to  the  axis  about  the  bond  mid  point  i s the  other  mode.  consequent  1 y , 
this  model  predicts  that  the  experimental  spectrum  of  the  solvation 
system  should  show  two  vibrational  peaks  associated  with  the  solvation 
structure.  \1 though  it  is  possible  that  the  second  peak  is  present 
(with  perhaps  ver\  weak  intensity),  the  authors  who  have  examined 
these  spectra  have  not  reported  am  additional  structure.  It 
seems  to  us,  t lie  re  fore,  that  in  spite  ot  the  suggestion  that  the 
solvation  number  ot  the  lithium  cation  in  I'MSO  is  two,"  this 
conclusion  does  not  soom  to  he  suhst ant i at ed  bv  the  far  infrared  spectrum. 

The  analysis  in  this  paper  is  based  on  the  assumption  that 
the  electrostatic  interact  ions  can  be  modelled  in  terms  of  the 
familiar  point  charge  interact  ions.  It  is  possible  to  test  the 
reasonableness  of  this  assumption  by  comparison  with  a potential 
energy  function  which  is  derived  on  the  basis  of  the  assumption  that 
the  charge  distributions  have  a diffuse  nature.  In  this  manner  it 
is  possible  to  enter  further  account  of  the  ionic  po 1 a r l : ab  i 1 i t v 
into  the  potential  energy  function.  It  will  be  evident  from  the 
form  of  tin'  potential  function  that  there  are  non  vanishing  electro 
static  contributions  to  the  force  constant.  Here,  in  the  remaining 
paragraphs  of  this  appendix,  we  consider  onlv  the  general  form 
of  the  interaction.  I he  derivation  of  the  electrostatic  cent r i but  ions 
to  the  force  constant  is  considered  later  in  appendix  4. 


ft 


The  form  of  the  interaction  between  a continuous  charge 

distribution  and  a point  dipole  located  at  the  centre  of  a conducting 

sphere  of  radius  a is  developed  with  the  use  of  the  Fourier 
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convolution  theorem  as  discussed  by  Si  1 verstone . “ The  general 
form  of  the  interaction  can  be  expressed  as 

V ( r i ? ) = -±-  jdJk  k'2p*oft  (k)PtUpole(k)expCik*ri2)  , (1.101 


where  ri?  is  the  vector  separation  between  an  element  of  charge 
in  the  continuous  distribution  and  the  point  dipole.  The  Fourier 
transforms  of  the  charge  densities  are  signified  by  the  p-quant i t ies . 
The  general  expression  for  the  transform  is 


p ( k 1 = d3r  p(r)exp(ik*r) 


(1.11) 


In  the  following  discussion  we  assume  that  the  continuous 

charge  distribution  is  modelled  as  a Slater-type  distribut ion, as 

has  been  discussed  in  connection  with  the  free  energy  of  ionic 

solvation  elsewhere  (cf.  ref.  25).  Thus,  the  transform  ol  the 
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charge  density  can  be  expressed  as"' 


P^offOp  - pV(p2  ♦ k 2 ) 2 


(1.12) 


corresponding  to  a coord i nate - space  distribution  of 
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(1.13) 


where 
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p = :/h 


(1.14) 


atnl  b is  t ho  effect  ivo  Bohr  radius  for  the  distribution.*' 

The  transform  for  the  dipolar  distribution  can  be  expressed 
a s 


d i po  1 


0(k)  = 4ni  lu/a).i  i (ka)  }'  YimU')Yimlk') 


(1.15) 


where  i i (ka)  is  the  spherical  Bessel  function  of  the  first  kind  of 

. .11 

o rde  r 1 . 

When  the  above  transform  expressions  are  substituted  into  the 
expression  for  the  interaction,  we  find 


V ( r , j ) = 


— j-liAiil — v (k  )exp  ( i k • r 
k7  (p7+k7)  7 ,m 


1 . 

(1.10) 


We  now  consider  the  expansion  of  this  function  with  respect  to 
two  arbitrary  centres  0 and  S (for  the  origin  of  the  centre  of 
svmmetrv  anil  the  location  of  the  solvent).  Thus,  if  R is  the  vector 
from  0 to  S and  r is  the  vector  from  0 to  the  location  of  the 
centre  of  the  charge  dens  i tv,  we  can  write  r. ? = R - r.  The 
exponential  function  in  (l.lo)  now  can  be  written  as  a product  of 
exponentials.  The  Ravleigh  expansion  can  be  used  for  each  of  these 
exponent i a 1 funct ions  : 1 1 
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As  a resul t , 
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Vho  angular  integrations  can  bo  carried  out  easily. 
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lhe  k space  radial  integral  lias  a pole  at  k = 0.  The  contribution 
ot  t h i s part  of  the  integral  to  the  whole  expression  is  identical 
to  the  ('arisen  Rushbrooke*  expression  as  wo  have  used  it  tor 
dipolar  systems.  Hie  remaining  terms  arise  from  the  pole  at 
k*ip.  \s  a result  , denoting  bv  \ ' all  terms  different  from  the 
t’a  r l son  - Rushbrooke  expansion,  we  can  write 
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1 he  angular  dependence  of  this  expression  is  the  same  as  the 
('arisen  Rushbrooke  expansion.  However,  there  is  an  important 
distinction  between  this  term  and  the  ('arisen  Rushbrooke  formula 
in  that  the  corotli  order  contribution  is  non  vanishing  in  its 
r dependence.  l'hus  , as  we  show  in  appendix  4,  there  is  an  electro 
static  contribution  to  the  force  constant. 
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Kith  l.  M.  McKinley,  Pepartment  of  Physics,  Oakland  University 


In  discussing  ionic  solvation,  so  tar  wo  have  concentrated 
on  1 v on  the  vibrational  degrees  of  freedom  which  exist  for  the 
system  of  the  ion  and  its  primary  solvation  shell.  This  subsystem 
i considered  to  be  located  inside  a cavitv  within  the  bulk  of 
the  solvent.  In  tins  section  we  examine  the  consequences  of 
adding  elect  rost.it  ic  image  effects  to  the  basic  cavitv  interaction 
t e rms . 

Iriedmaii  ' has  considered  the  image  approximation  to  the 
dielectric  reaction  field.  His  approximation  methods  are  particularly 
useful  for  the  cons i derat  ion  of  complicated  collections  of  solvent 
and  ions  within  a solvation  cavitv.  Although  Iriedmaii  considered 
the  image  approximat ion  to  the  monopole  t icld  in  detai 1 , he  onlv 
quotes  the  results  for  the  dipolar  field.  lor  the  sake  ot  complete- 
ness in  our  discussion,  we  outline  the  electrostatics  ot  the 
dipolar  reaction  field  and  the  application  of  the  image  approximation 
to  it. 


Consider  a single  charge  q located  at  s on  the  c axis  inside 
a spherical  cavitv  of  radius  a.  The  potential  due  to  this  charge 
inside  the  cavitv  is  given  by-  • " 
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where  P lx)  is  the  1 egemlre  polynomial  or  order  n. 
spherical  cavitv,  the  potential  is 


Out  s i vie  t he 


out,r,0)  = ^TTrnrrrr  pfrnr  {Vcos0)- 


(2.:) 


ho  now  imaSine  ■'  single  image  charge  q-m  placed  at 
lhe  potential  produced  hv  this  charge  at  r<a  is  given  by 
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would  like  to  make  the  identificati 
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!his  i dent i t icat ion  will  be  valid  for  all  sufficiently  large  values 
01  f*  *hus,  r.ffl  and  q.^  can  be  identified  as 


r . * a * / s 
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(2.5) 


and 
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A 1 1 e i na t i ve  1 v , 1-  r iednian’’ ^ lias  noted  that 
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Thus,  the  identifications 


(2.5)  and  (2.o)  amount  to  the 


as sump t i on 
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that  only  the  k = 0 term  in  (2.7)  needs  to  he  considered.  It  is 

equivalent  to  say  that  for  sufficiently  large  values  of  C,  the 

k = 0 term  is  the  dominant  contribution.  It  should  be  noted,  however, 

that  also  for  sufficiently  large  values  of  the  dielectric  constant 
t e rm 

c the  k=0/also  suffices  regardless  of  the  size  of  C. 

If  we  consider  only  the  k = 0 contribution  to  (2.7),  eqn  (2.3) 
can  be  written  as 


*im  = - —pr  ( a / s ) ^ — J(r/a:s)erc(cos 0)  (2.8) 

im  1 1 (a?/s)  C C 

which,  as  Friedman  has  pointed  out,"  is  equivalent  to  the  Laplace 
expansion  of 


‘lin/frTiJ 

for  r < r . m . 

In  order  to  consider  a dipole  located  in  the  cavity  and  its 
image,  we  next  consider  the  superposition  to  two  solutions: 
one,  0 1 is  the  potential  produced  by  q j = +q  located  at  : = s,  and 
two,  is  the  potential  produced  by  q2  = -q  located  at  ; = s-<5. 

The  potential  inside  the  cavity  produced  by  q»  is  given  by  eqn  (2.1). 
The  potential  produced  by  q2  is 
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The  total  potential  inside  the  cavity  is 
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\s  before  with  eqn  (2.31  ami  the  identification  (2.4),  here  we 
want  to  establish  a similar  means  of  identifying  the  dipole  image 
and  location  in  terms  of  the  dipole  moment  and  location  inside  the 
cavitv.  Ihus,  in  eqn  (2.13)  we  add  and  subtract  the  same  quantitv: 
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The  second  term  is  identified  as  the  image  dipole 
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As  before  for  the  monopole  case,  in  the  limit  of  large  C or  e,  the 
location  of  the  image  is  given  bv 


i m 


(2.19a) 


and  the  dipole  moment  of  the  image  is  given  bv 


q.  3 = (a /si  \iS  . 
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finally,  the  third  term  in  12.1")  is  identified  as  an  added  image 
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Tho  purpose  of  the  imago  approximation  is  to  simplify  the 
ca  leu  1 at  ions  involving  the  use  oi  the  reaction  field.  lhis  is  done 
hv  replacing  the  complicated  solutions  to  the  I.aplace  equation 
h\  simpler  expressions  which  involve  the  interaction  of  discrete 
chare.es  or  multipoles.  The  image  approximation  amounts,  simply, 
to  the  assumption  of  the  k=0  term  in  eqn  (2.”’). 

l'he  energy  of  interaction  of  a charge  with  its  own  image,  or 
more  exact lv  its  reaction  field,  defines  the  self  energy  ot  that 
charge  with  respect  to  the  dielectric.  The  sell  energy  is 
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l'he  interaction  of  any  other  charge  with  this  image  charge  can  be 
expressed  as 


Thus,  we  are  now  in  a position  to  calculate  the  image  contribution 
to  the  ionic  potential  energy  function  for  the  solvated  ion. 

In  this  manner  we  can  approximate  the  contribution  ot  the  remainder 
of  the  solution  svstom  to  the  vibrational  modes  of  the  ion  and 
its  primarv  solvation  svstem. 

In  the  following,  we  illustrate  the  image  approximation 
with  the  use  of  the  tetrahedral  solvation  model. 
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Within  the  solvation  cavity,  the  terms  which  contribute  to  the 
breathing  mode  force  constant  are  the  following:  (1)  the  ion-dipole 
interaction,  (2)  the  dipole-dipole  interaction,  (3)  the  short- 
ranged  interactions  (repulsions)  which  operate  between  the  species 
in  a pair-wise  sense,  and  (4)  the  interaction  between  the  solvent 
molecules  and  the  cavity  surface.  Outside  the  cavity,  in  the  bulk 
of  the  dielectric,  the  various  image  quantities  are  located. 

We  consider  first  the  electrostatic  interactions  which  operate  between 
species  within  the  cavity. 

The  ion-dipole  interaction  has  been  discussed  earlier.  For  the 
tetrahedral  solvation  it  is  simply 


Vid(c) 


4ne/R2 . 


(2.22) 


We  next  calculate  the  dipole-dipole  interaction  which  operates 
between  each  pair  of  solvent  molecular  dipoles  within  the  cavity. 

To  carry  out  this  calculation,  we  assume  that  the  dipoles  can 
he  regarded  as  ideal  point  dipoles.  Thus,  the  general  form  of  the 
interaction  is  given  by 


Vjd  = Mi  * M 2 /R 3 - 3 (m i • R) (m 2 * R) /R5  (2.23) 

and  the  vector  R is  given  by 

R = R i - R 2 . (2.24) 

The  vectors  Ri  and  R2  are  referred  to  a common  origin  at  the  centre 
of  symmetry  for  the  system. 
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' wa  1 1 1 + exp  [ (ft-A)  fix  1 

where  A is  the  cavity  radius  and  a is  the  half  width  of  the  barrier. 

The  thickness  of  the  harrier,  and  the  extent  of  its  influence  into 
the  cavity,  is  determined  hv  the  magnitude  of  the  constant  a.  Tor 
a value  of  V0  of  the  order  of  several  cV  (a  magnitude  which  corresponds 
roughlv  to  the  solvation  well  depth),  the  effect  of  the  potential 
term  (2. do)  is  very  small.  Its  effect  on  the  breathing  mode  force 
constant  is  even  less. 

The  contribution  of  these  cavity  potential  terms  to  the 
solvation  shell  harmonic  breathing  mode  force  constant  is  easily 
found  by  direct  differentiation.  Thus, 

k = -24  iiil  + 45/571  +20  R/p 2 +4^-  cxPl  ( R o ~.A  1 /aj_~_l (2.27) 

c Ro  Ro  a2  (exp [ (R0-A)/a]+l) 3 

where  R and  p,  associated  with  the  repulsion  interaction,  are  defined 
in  the  text:  cf.  eqn  (3). 

We  now  turn  to  a consideration  of  the  image  contributions  to 
the  potential  energy  function  and  to  the  force  constant.  These 
terms  show  a complicated  structure,  and  their  ultimate  effect  on 
the  force  constant  is  small.  Thus,  we  will  show  that  they  generally 
need  not  be  considered  further,  especially  when  one  is  interested 
only  in  the  vibrational  structure  of  the  system. 

There  are  a number  of  contributions  to  the  potential  energy 
function  for  the  system  which  arise  from  the  interaction  of  the  ion 
and  the  solvent  dipoles  with  images  in  the  bulk  of  t lie  solution. 

Two  of  the  interactions  yield  no  contribution  to  the  harmonic  force 


■I  0 

force  con  s t ;iii  t , although  t hev  in  a y In'  of  considerable  importance* 
in  matters  concerned  with  solvation  a ml  free  energy  changes. 

In  particular,  the  ionic  self  energy  (the  interaction  of  the  ion 
with  its  i matte)  tlepemls  on  1 v upon  the  solvation  shell  ratlins  A. 
•Mso,  the  interaction  between  the  ion  insitle  the  cavitv  ami  the 
images  of  the  tlipoles  yields  no  contribution  to  the  force  constant. 
These  terms,  ami  others,  will  he  illustrated  in  the  following 
pa  rag raph s . 

For  the  classical,  non  spat  ini lv  ami  temporal Iv  disperse 
dielectric,  the  ionic  self  energy  is  given  by  the  Born  expression 


is) 


(1  l/. r,x 


U.  :s) 


n t lu*  image  approx imat ion  of  this  expression,  we  t i ml 


J(* 


Y 


i is) 


1 r 1 e- 

? r*i  \ • 


i : . : ;>) 


For  DMSO,  for  example,  the  image  approx i mat i on  differs  from  the 
Born  self  energy  hv  about  For  water,  the  image  approximation 

differs  from  the  Born  self  energy  by  only  1.31.  Thus,  for  polar 
solvent  systems,  the  image  approx imat ion  is  reasonable.  In  any 
vase,  it  is  clear  from  the  form  of  oqn  (J.JS)  and  U’.-’F)  that  there 
can  be  no  contribution  arising  from  these  terms  to  the  force 
con  s t ant . 

In  the  image  approximation,  the  sell  energy  ol  the  dipole  is 
relatively  easy  to  determine.  Making  use  of  the  location  of  the 
dipole  image,  we  can  write  the  vector  separation  between  any  dipole 
and  its  image  as 
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Rw  = (1  - A7/R7)R  (2.30) 

where  R is  the  vector  drawn  from  the  centre  of  symmetry  to  the 
dipole  in  the  cavity.  The  image  dipole,  and  its  source,  lie  along 
the  vector  R (and  its  extent  ions)  , but  these  dipoles  are  directed 
in  the  opposite  sense.  from  the  general  form  of  the  d ipole -d ipole 
interaction,  we  find 


Vdd(i) (s) 


H ? A 1 

(A7-R7) 3 


(2.31) 


As  we  discussed  above,  and  as  Friedman  also  notes,  co-located 
with  the  dipole  image  is  an  effective  charge  with  a magnitude  given 
by 


(2.32) 


The  interaction  between  the  source  dipole  and  this  part  of  the  image 
can  be  written  as 


r - 1 u 7 


d c ( i ) ( s ) r+1  ,A,  (1-rVA7) 


(2.33) 


When  (2.33)  is  combined  with  (2.31),  we  can  write 


V 


dd(s) 


l » - 1 n?  l + x7 

'7  7 TT 


a 3 ( l - x 7 ) 3 


(2.34) 


in  which  x is  defined  bv 
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x = R/A . 


(2.35) 


This  quantity  is  always  loss  than  or  at  worst  equal  to  unity. 

By  means  of  a direct  differentiation  (with  respect  to  R) 
of  (2.34)  wo  find  that  this  term  contributes 


L . . ' ^_4  111  _L2 (’Ox*  * 1 1 x ? - 1 1 (2.30 

to  the  force  constant.  I ntorost i ng 1 y enough,  this  term  contributes 
maximally  (viz.,  zero)  to  the  force  constant  when  R = 0.28A. 
Conorallv,  however,  a solvation  radius  as  small  as  0.28A  is  unlikely. 

The  interaction  between  the  ion  and  the  dipole  images  is 
easily  found.  It  is  given  by 


V i d ( i ) 


(2.37) 


As  remarked  earlier,  this  term  does  not  make  a contribution  to  the 
force  constant  (it  is  independent  of  R) . 

Next,  we  consider  the  interaction  between  the  ion  and  the 
effective  charges  associated  with  the  dipole  images.  This  term 
is  given  simply  by 


,,  i-l  ou 

vi  i m = rrr  att- 


(2.38) 


The  interaction  between  a dipole  and  the  dipolar  image  which 
arises  from  another  dipole  in  the  cavity  is  fairly  complicated. 


The  derivation  is  sketched  as  follows.  The  distance  between  any 
dipole  and  the  image  of  any  other  dipole  is  given  by 


R i 2 = [3(A‘*  + R4)+2A2R2]  1^2  / /Jr  . 


(2.39) 


for  a system  of  tetrahedral  solvation,  the  scalar  product  of  the 
dipole  with  the  ot her -d ipo le - image  is 

M*U2  = fr}(A/R)3  4 • (2-40) 

The  scalar  product  of  the  image  dipole  moment  with  the  radius  vector 
between  the  dipole  and  the  other  image  is 

uI-Ri2  = f^j-yR(A2/R2*l/3)  (A/R)3.  (2.41) 


Similarly,  for  the  cavity  dipole  we  find 


y 2 * R 


uR 


(A/R)2+l 


(2.42) 


The  general  d ipol e -d i pole  interaction  is  now  evaluated  to  give 


Vdd ' ( i ) 


2/3 


c - 1 

c + l 


y 2 A 


3Ai*  + 14A2R2  + 3R4 
(3A4  + 2A2R2  + 3R‘,)5/2  * 


(2.43) 


Finally,  in  a similar  manner,  we  find  that  the 
a dipole  inside  the  cavity  and  the  effective  charge 
the  dipolar  image  arising  from  a separate  dipole  in 
given  by 


interaction  between 
co- located  with 
the  cavity  is 
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4 4 - 


do  l i 4 
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rn 


x 7 ( x ? ♦ 3 1 


■\*  (3  + 2x’  + 3x>'771 


(2.444 


The  quant  i tv  \ is  defined  hv  oqn  (2.354. 

I rom  all  those  cent  r i hut  i ons  to  tho  potential  energy  l unction, 
wo  t i ml  tho  imago  contrihution  to  tho  force  constant  is  given  by 
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( 2 . 4 ' 4 


The  quant i tv  \0  is  R0/A  whore  R0  is  the  equilibrium  solvation  radius. 
I'ho  difference  between  the  cavitv  contributions  to  the  force 


- - ■ ^ - - - - ' 
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constant  ami  tlu*  total  form  of  the  force  constant  as  a function 
of  the  solvation  rail i us  of  a tetrahedral  solvation  system  is 
illustrated  in  figure  1.  The  difference  is  not  significant. 

The  image  contributions  to  the  force  constant  for  t ho  breathing 
mode  are  small,  and  probably  can  be  ignored  in  most  cases. 


Appendix  3:  The  l.aplace  expansion  of  the  exponential  repulsion 


As  we  indicated  earlier  in  this  paper,  in  order  to  be  able 

properly  to  consider  the  vibrational  modes  associated  with  the 

motion  of  an  ion  within  its  solvation  shell,  it  is  necessary 

to  make  use  of  an  appropriate  form  of  the  potential  energy  function. 

8 

for  the  electrostatic  terms  the  Carlson -Rushbrooke  expansion  was 
used.  The  scalar  exponential  repulsion  which  operates  between 
the  ion  and  any  of  the  solvent  molecules  can  be  expanded  in  a 
similar  appropriate  manner.  Tt  is  the  purpose  of  this  appendix 
to  outline  the  construction  of  this  expansion. 

The  scalar  exponential,  Born  repulsion  has  the  form 

V(ri2)  = B exp  [ - (r j 2-a)/p]  . (3.1) 

The  expansion  of  this  function  in  terms  of  the  coordinates  r>  and 

r<  (where  r>  is  the  greater  of  ri  and  r2)  is  carried  out  in  a 
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manner  similar  to  that  of  Ruedenberg,  and  others,4-  ’ for  the 
Laplace  expansion  of  the  free-space  Green  function  (viz.,  l/ri2). 

Iv'e  make  use  of  the  Fourier  transform  of  eqn  (3.1)  which  is 
given  by 

v(k)  = 8nR  exp(a/p) . (3.2) 

p(l/p2+k2)2 

G i ven 

f(k)  = jd3r  F (r)exp [ i k • r ] , (3.3) 

it  is  well  known  that  one  can  write 
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(3.4) 


F(r-R) 


1 

C2n ) 3 


d3k  f (k) 


exp [ - ik • (r -R) ] 


Thus,  for  example,  if  f(k)  = 4n/k2  and  the  exponential  function  in 
eqn  (3.4)  is  written  as  a product  of  two  exponential  functions, 
then  the  Rayleigh  expansion  can  be  used  in  place  of  each  exponential. 
When,  in  the  resulting  expression,  the  various  integrations  are 
carried  out,  one  recovers  the  familiar  Laplace  expansion  of 
1/ I r-R | . 

We  now  make  use  of  v(k)  in  place  of  f(k)  in  eqn  (3.4)  to 

write 


V(r-R) 


16  5.  axpta/p)  l vj  Cr)Yt  (R) 

k ,m 


3e  (kr)j»  (kR) 

dkk2  — 

(1/ p2+k2) 2 


(3.5) 


There  is  no  pole  at  k=0.  Therefore,  we  need  only  evaluate  the 
residue  of  the  radial  integral  at  k=i/p.  The  general  form  of 
the  expansion  of  (3.1)  now  is  written  as 


V (r-R)  =-4iiRexp  (a/p)  l Y£m 

£ ,m 

+ (r</p)iJl  + 1 (r</p)k4(r>/p) 

where  r>  is  the  larger  of  r and 
The  expansion  of  the  Morse 
eqn  (3.6).  The  Morse  potential 


/v 


(r)Y*m(R)  (2a+l)il(r</p)k1(r>/p) 


- (r>/p)ia(r</p)kjl+1(r>/p) 


(3.6) 


R. 

potential  is  merely  an  application  of 
function  is  given  by 


I'{  1 - exp  ( - ( r - ) /a]  ) 2 


= D { 1 - ’ exp [ - (r - r^) /a  1 + exp [ - J ( r - ) /a ] ) . (3." 

The  Fourier  transform  of  the  first  term,  H,  yields  the  Dirac  delta 
function.  Consequently,  the  inverse  transform  yields  the  constant 
p in  return.  The  remaining  terms  in  eqn  (3.  ""I  behave  in  a manner 
identical  to  the  Born-type  repulsion.  Thus,  the  expansion  of  these 
terms  is  accomplished  directly  with  the  use  of  eqn  (3.M.  When 
the  contributions  from  the  two  exponential  terms  are  grouped,  the 
form  of  the  expansion  given  in  eqn  ("")  results. 
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Appeiuiix  4:  Harmonic  force  constants 

(liven  a potential  enerjjv  function  for  a system,  the  second  order 
term  in  the  l’aylor  series  expansion  defines  the  harmonic  oscillator 
potential  energy  function.  I'h  i s second  order  term  is  determined 
hv  the  d i f ferent i a t i on  of  the  original  function,  followed  by  the 
evaluation  of  the  resulting  expression  in  terms  of  the  equilibrium 
values  of  the  coordinates.  This  is  t lie  obvious  route  to  the  force 
constant  quantities. 

However,  for  a number  of  complicated  functional  forms  of  the 
potential,  the  use  of  differentiation  can  be  cumbersome.  In 
particular,  when  one  seeks  the  form  of  the  oscillator  potential 
which  operates  on  the  ion  inside  the  rigid  cage,  the  differentiation 
of  the  complicated  angular  dependencies  is  not  necessarily 
straightforward.  Indeed,  if,  instead  of  the  symmetric,  regular 
systems  considered  in  this  paper,  one  were  to  consider  systems 
in  which  the  solvent  molecules  were  arbitrarily  placed  in  the  solva- 
tion cavity,  the  determination  of  the  oscillator  potential  would  be 
extremely  difficult.  Such  differentiations,  as  would  be  carried 
out  by  this  method,  could  be  pursued  with  the  use  of  the  formula 
derived  by  Rose.'* 

On  the  other  hand,  as  we  have  shown  in  Appendices  1 and  3, 
often  the  potential  energy  function  can  be  expressed  in  terms  of 
the  Fourier  convolution  integral.  In  this  case  it  is  simpler  to 
differentiate  the  integral.  The  analysis  then  proceeds  with  the 
evaluation  of  the  integral.  IV e will  show  that  in  this  manner  the 
calculation  often  can  be  simplified. 
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In  the  following  paragraphs  we  will  derive  the  integral  repre- 
sentation o f the  second  order  Taylor  series  term.  The  first  order 
term  can  he  found  in  a similar  manner.  However,  here  only  the 
resu It  will  he  quot ed . 

Given  a function  which  can  he  expressed  in  the  form 


G(r) 


1 ( r 1V  \ 


(rl  , 


14. n 


we  wish  to  find  an  integral  representation  of  the  seconw  order 
term  which  can  he  expressed  as 


t ? ( r ) = 1 i m 4 ( r • V ) 2 G (a  1 , (.  4 . J ) 

— w - » l 

rro 

where  r0  is  the  equilihrium  value  of  r.  The  notation  implies  simple 
that  in  G(.r),  r is  replaced  hv  a,  but  the  functional  form  remains 
unchanged.  The  differential  operators  now  operate  only  on  the 
function  G.  As  a result,  it  is  possible  to  write 


t : i r 1 


1 im  i rr :V  V G (a) 

..  » » kl  J 

a 'To 


1.4 . 31 


As  the  effect  of  the  differential  operators  now  has  been  established, 
we  can  drop  the  distinction  between  the  coordinates  a and  r without 
risk  of  confusion. 

The  l-'ourier  transform  of  G(rl  is  written  as  f(k)Y.  (k)  . Thus, 
equn  (.4.3)  can  he  written  as 
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Kith  the  interaction  between  the  continuous  charge  distribution  and 
the  point  dipole  which  was  discussed  in  Appendix  1. 

The  I'ourier  transform  of  the  exponential  repulsion  is  given 
In-  eqn  (A. 21.  We  are  interested  in  the  harmonic  term  witli  respect 
to  t he  r variable  in  the  general  expression  v(r-R):  cl'.,  eqn  (3.4) 
and  (3.S).  1'hus,  we  write 


V t r - R ) 


1 

(Jn)  ' 


d ' k | v ( k ) e x p ( i k • R ) ] e x p ( - i k • r ) . 


(4  . 1 2) 


P i f f e ren t i a t i ng  this  expression  twice  with  respect  to  r yields 


V V V ( r R) 
r r 


( 2 n ) 3 


d'k  kk  [ v (k ) exp ( i k • R) | exp ( - i k • r) • (4.13) 


The  integral  in  (4.13)  has  the  same  form  as  the  integral  expression 
in  eqn  (4.4).  Subsequent  evaluation  of  eqn  (4.13)  proceeds  with 
the  expansion  of  the  exponential  term.  Thus, 


| v ( k ) exp  ( i k • R ) | = )"  f^(k)YX|(k) 

*,U  11 

32 n •’  Rexp  (a/ p) — — )'  iXY  (R)Y’  (k) 

P ( 1 / p 7 + k 7 ) 2 \.n  M U 

\ix(kR).  (4.14) 


The  lj  Iro)  defined  by  eqn  (4.8)  now  has  the  specific  form 
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n 


■ilfo)  = exp  (a/p)  [ iX)\  (R) 


A ,u 


Am 


a’  .i , (kr0)  j . (kR) 

dk  k4  — — - 

o (l/p2+k2)2 


(4.151 


The  integral  is  evaluated  for  r„  < R,  which  will  always  he  true 
(the  limit  is  taken  as  rovO).  There  is  no  pole  at  k = 0 unless 
A''l>4,  which  will  not  happen.  As  a result,  we  find 

[°°  j,  (kr„l  j (kRl  , 

K‘  (tvTiee  * 1 * 

+ (r/o)  i j + j (r/p')kx  (R/p)  (R/p)  i t (r/p)kA+  j (R/p)  } . (4. Ip) 

1 lie  equilibrium  value  of  r0  is  0.  As  a result,  the  integral  (4.  Ip) 
vanishes  for  all  values  of  1.  } 0.  Consequently,  there  is  only 
one  term  1 0 ( r0 ) : 

lo(ro)  = iS  it 2 ( M/  P 2 ) exp  ( a/  p ) [ V x ( R)  { ( m\  ) k .(R/p)  - (R/p)k.  .(R/p)} 

\ , n ' + 1 

(4.17) 

However,  the  first  term  in  (4.11)  is  I x ( r 0 ) . That  is,  the  first 
term  in  (4.11)  corresponds  to  l.  = .\.  hut,  since  1=0,  only  one  term 
* 11  (^-1  ) contributes.  Thus,  the  first  term  in  (4.11),  evaluated 
for  this  system,  is 


t?(r)  = r 2 ( R/b  p 2 ) exp [ - ( R - a ) / p ] ( 1 - 2 p/R ) . 

hhen  this  quantity  is  evaluated  at  the  svstem  equilibrium  point 
of  R = a = R , we  find 

t?(r)  * r 2 (R/pp 2 ) (1  - 2p/R  ). 

e<i 


(4.10) 
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The  second  term  in  eqn  (4.11),  evaluated  for  this  repulsion, 
has  the  form 

1 2 (r)  = — ( R/ p 2 )exp (a/ p)  l /2TTT(A200 | 00) (A2u-U | 00) 

5/5  A , u 

*Y*^  (R)  Y*  (r)  {(R/p)kA  + 1 (R/p)  - (3  + A)kx(R/p)}  . (4.20) 

Clearly,  the  sum  collapses  to  one  term,  the  term  for  A=2.  For 
any  system  with  simple  cubic  symmetry,  the  summation  of  the 
repulsion  terms,  in  the  same  sense  as  the  summation  of  the  dipole 
terms  in  the  electrostatic  case,  yields  a zero  contribution. 

Hence,  only  the  zeroth  order  term  survives;  that  term  is  given  by 
eqn  (4.19). 

The  application  of  this  analysis  to  the  Morse  potential  is 
striaght forward . We  proceed  from  eqn  (4.19). 

The  Morse  potential  can  be  seen  to  yield  two  terms  to  the 
force  constant.  The  force  constant  contribution  associated  with 
the  attraction  term  we  label  k+,  the  term  associated  with  the 
repuslion  contribution  we  label  k_.  Thus, 

k + = - l!^i(i-2/aR)  (4.21) 


and 


k 


4Da2 


(1  -1/aR)  . 


(4.22) 


The  sum  of  these  two  terms  gives 
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k = 2Daz/3  (4.23) 

We  turn  now  back  to  a consideration  of  the  interaction  between 
the  soft  ionic  charge  distribution  and  the  collection  of  point 
solvent  dipoles  which  was  discussed  in  the  first  appendix.  At 
this  point  we  need  to  assess  the  contribution  this  type  of  electro- 
static interaction  makes  to  the  force  constant. 

In  eqn  (1.1b)  we  expand  the  exponential,  exp(ik*ri2)»  in  terms 
of  r and  R:  ru  = R-r.  Therefore,  consistent  with  eqn  (4.4),  we 
write 


j i (ka) 


m k2(p2+k2)2 


exp ( i k • R) 


■ "■4-vi-iif>  i i i (iwBiMj)1/2(ijo(i|>ouijmt|>i,H>Yii.(°> 

X 1 1 m i k v y ' 


X , p m i , k 

j i (ka) j . (kR) 


xY..(R)Y,  (k)  * 

jk  An1  1,2  .1,2-v  2 


kz(p2+k2) 


(4.24) 


l f,  (MY  Ck) 
A ,y 


The  radial  quantity  Ij  ( r o ) is 


IjCr.)  - ^4^  l l 

m j ,k 


[ 3(21+1) 


4 7i  ( 2 A + 1 T 


1/2 


(1 jOO | AO) (1 jmk | Ay)  i 


; j 


*Ylm^VR) 


r j i (ka)  j . (kR).i . (kT0) 

1 dk  k2 i t 


(P2+k2)2 


(4.25) 
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IVe  now  evaluate  the  first  term  in  eqn  f4.ll).  Here,  as  with 
the  repulsion  term  examined  above,  we  note  that  I,  (r0)  vanishes 
as  ro  * 0 for  I.  + 0.  Note,  however,  that  the  integral  in  eqn  (4.25) 
does  not  depend  upon  \.  Nevertheless,  the  first  term  in  eqn  (4.11) 
requires  that  L=A=0.  IVe  find  that  the  first  term  in  (4.11) 
yields 


t ? ( r ) = 2 tt  r 2 


.>a 


IV, , (a ) Y 


1 -1 


(R) 


' 1 o 1 o ^ 


Yl-1  (a)Yn  (R)] 


x { 3i , (ap)k i (Rp)  + (ap) i 2 (ap) k i (Rp) 


For  a tetrahedral  svstem,  this  becomes 


(Rp) i i (ap) k 2 (Rp) ) . 

(4.20) 


t|Cr) 


( 3 i | (ap)k  j (Rp)  + 


(ap) i 2 ( a p ) k i (Rp) 


- (Rp) i i (ap)k2 (Rp) ) . 


(4.27) 


There  is  another  possible  contribution  to  the  force  constant 
for  the  case  when  l.  = 0,  ? = A = 2.  However,  it  is  relatively  easily 
shown  that  the  form  i<  similar  to  the  Carl  son - Rushhrooke  expression 
for  C=2.  Thus,  for  simple  cubic  systems,  these  additional  terms 
vanish  in-  s v mm  e try. 

The  behaviour  of  the  force  constant  contribution  (4.27)  as 
a function  of  solvation  radius  is  shown  in  Table  10.  The  first 
three  values,  which  correspond  roughly  to  the  solvation  radii  for 
lithium,  sodium  and  potassium,  can  be  seen  to  be  decreasing. 

On  the  other  hand,  the  last  two  values,  which  correspond  approximately 
to  the  solvation  radii  for  rubidium  and  cesium,  increase.  The 
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contribution  to  the  force  constant  has  a negative  value  as  a result 
of  the  attraction  between  the  soft  charge  and  the  point  dipoles. 

The  repulsion  contribution  will  continue  to  be  a positive  valued 
contribution  to  the  force  constant.  This  term  generally  will 
decrease  with  increasing  solvation  radius.  As  a result,  the  tirst 
three  force  constants  for  lithium,  sodium,  and  potassium  will  be 
expected  to  decrease,  as  observed.  However,  because  the  electro- 
static contribution  increases  from  rubidium  to  cesium,  the  overall 
force  constant  for  these  two  species  should  at  worst  decrease 
slower  than  is  the  case  for  lithium  to  potassium,  and  at  best,  they 
should  increase,  as  t lie  experimental  data  indicate. 

This  soft  charge  contribution  to  the  force  constant  does 
support  the  results  of  our  analysis.  The  softness  oi  the  charge 
distribution  is  connected  to  the  polarizability  of  the  atom  or  ion. 
Hence,  the  result  we  find  for  this  e lec t r i st at i c contribution 
i s not  su rp r i s i ng . 

finally,  to  complete  this  discussion  of  force  constant  problems, 
we  quote  the  integral  form  of  the  first  order  Tailor  series  term. 

I'h  i s term  generally  is  needed  any  time  a particular  distribution 
of  solvent  dipoles  deviates  from  the  equilibrium  symmetry  and 
arrangement  of  some  reference  state. 

The  first  order  contribution  to  the  Tavlor  series  expansion 
is  wr i 1 1 on  as 
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fable  3:  Solvation  radii 
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Table  4_:  Calculated  force  constants  according  to  eqn  (9a)  and  (9b). 
The  results  are  presented  in  the  following  manner:  the  first  number 
is  the  force  constant;  the  first  number  in  pa  rent  he sc s is  the 
calculated  vibrational  frequency  (in  cm  *)  using  the  reduced  mass 
associated  with  the  coordination;  the  second  number  in  parentheses 
is  the  calculated  frequency  assuming  a rigid  solvat  ion  cage  (c  = >»)  . 
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